Establishment of a zebrafish inbred strain, M-AB, capable of regular breeding and genetic manipulation

Inbred strains of organisms are genetically highly uniform and thus useful for life science research. We have previously reported the ongoing generation of the zebrafish IM strain from the India (IND) strain through full sib-pair mating for 16 generations. However, the IM fish laid a small number of offspring and had a short lifespan, implying the need for discreet care in breeding. Here, we report the subsequent establishment of IM strain as well as the generation of a new inbred zebrafish strain, Mishima-AB (M-AB). M-AB was derived from the *AB strain by full sib-pair mating for over 20 generations, which fulfills the general criterion for the establishment of an inbred strain. In contrast to the IM case, maintenance of the M-AB strain by sib-pair mating required almost no special handling. Genome sequencing of IM individuals from the 47th generation and M-AB individuals from the 27th generation revealed that SNP-based genomic heterogeneity across whole-genome nucleotides was 0.008% and 0.011%, respectively. These percentages were much lower than those of the parental IND (0.197%) and *AB (0.086%) strains. These results indicate that the genomes of these inbred strains were highly homogenous. We also demonstrated the successful microinjection of antisense morpholinos, CRISPR/Cas9, and foreign genes into M-AB embryos at the 1-cell stage. Overall, we report the establishment of a zebrafish inbred strain, M-AB, which is capable of regular breeding and genetic manipulation. This strain will be useful for the analysis of genetically susceptible phenotypes such as behaviors, microbiome features and drug susceptibility.

the establishment of a separate inbred strain that can be maintained by regular care has been a desire of the research community.
In this study, we attempted to generate a new inbred strain that is better at breeding than the IM strain.To eliminate unwanted genetic variants that cause inbreeding depression, we employed the *AB strain, which was generated after efforts were made to remove harmful variants through the selection of good homozygous gynogenetic females that produce haploid embryos showing normal early embryogenesis 12 .We carried out full sib-pair mating for more than 20 generations and established an inbred strain named Mishima-AB (M-AB).Genome sequencing of IM and M-AB individuals revealed that the genetic heterogeneities of these strains were low.Embryos of the M-AB strain fish were robust in development, and applicable to microinjections at the 1-cell stage.Collectively, the inbred M-AB strain is amenable for breeding and useful for genetic studies in zebrafish.

Subsequent establishment of the inbred IM strain
The IM strain was generated by repetitive sib-pair mating of the IND strain for 16 generations 11 .We continued sib-pair mating for a total of 20 generations and here declare the establishment of the IM strain as an inbred strain (Fig. S1A).Further sib-pair mating of the IM strain fish was performed, reaching the 50th generation in 2023 (Fig. S1B-D).During this sib-pair mating, inbreeding depression was observed in many branches.The embryonic lethal phenotype with a small head, an underdeveloped lower jaw and poor circulation, which we have previously referred to as the Tübingen-Mishima (TM) phenotype 11 , was observed in the clutches of the I 16 -3, I 20 -5, I 22 -9, I 23 -4, I 23 -5 and I 23 -10 pairs (Table S1).The total number of eggs per clutch in generation 20 was less than half of that in generation 0 (Fig. S2A).The efficiency of fertilization (number of fertilized eggs/number of total eggs) declined to 50% in the 19th generation (Fig. S2B,C).The survival rate of offspring at 3 days post fertilization (dpf) also declined after the 14th generation (Fig. S2D,E).Taking into account these inbreeding depression data, we generated various branches of the IM subfamily after the 15th generation and selected a few good branches for subsequent maintenance.This labor-intensive inbreeding improved the number of eggs, efficiency of fertilization and survival of offspring from the 23rd generation.The 50th generation of the IM strain is currently maintained but has a low fertility rate and a small number of embryos (Table S2).

Generation of the inbred M-AB strain
To attempt to generate a new inbred line with lower inbreeding depression, we employed the *AB strain, which was generated from the AB strain by selecting good parthenogenesis-derived females 12 .A course of sib-pair mating started in 2010 followed almost the same procedure as that used for the IM strain except for maintaining a smaller number of branches (Fig. S3A).During consecutive inbreeding, we neither observed the fragile TM phenotypes nor found a reduction in egg number or fertility rate.However, the *AB-derived strain fish often gave rise to male-biased progeny that resulted in the discontinuation of several branches, i.e., AB 19 -7 and AB 19 -10 pairs (Fig. S3B).This tendency of male bias was mitigated by switching the larval diet from paramecia to rotifers and by lowering the growth density of larvae and juveniles in a tank (see "Materials and methods").Sib-pair mating was continued for over 20 generations, and we declare the establishment of the M-AB strain as an inbred strain.Consecutive sib-pair mating of the M-AB strain reached the 30th generation in 2023 (Fig. 1A).
The M-AB fish looked healthy (Fig. 1B) and were maintained by regular fish care but with a concern of a male-biased tendency that appeared in several generations.The number of fertilized eggs (~ 200 embryos), the www.nature.com/scientificreports/efficiency of fertilization (~ 90%) and the length of the body at 30 dpf were almost comparable between M-AB in the 27th generation and the parental *AB strain (Fig. S4A-D).

Phylogenetic analysis of the M-AB and IM strains
To assess the phylogenetic relationship of the inbred M-AB and IM strains with their parental strains, we performed whole-genome sequencing (WGS) of three individual fish of the 27th generation of M-AB (AB 27 -4), the 47th generation of IM (I 47 -2) and our lab stocks of *AB and IND.We also obtained genomic information for AB (AB_SRA) and TU (TU_SRA) from the NCBI Sequence Read Archive (SRA) as well as the zebrafish genome assembly of the TU strain (TU_GRCz11) from the Genome Reference Consortium 13 .We used TU_GRCz11 as a genome reference to identify SNPs in our 12 fish and the 2 archived sequences (Table 1).The total number of SNP positions found in the 14 total genomes was 24,312,666 (18.1% of the zebrafish genome).As expected, the number of SNPs found in the TU_SRA genome was lower than that in the other genome.Based on these SNPs, we created a phylogenetic tree of the 14 genomes by the maximum likelihood method using the standard transversional model 14,15 .The 3 individual fish of each strain were highly similar in each subgroup (Fig. 2).The M-AB and IM subgroups showed the closest relationship with their parental *AB and IND strains, respectively.These data indicate that M-AB and IM were phylogenetically isolated from their parentals by repetitive sib-pair mating, recapitulating the establishment of the inbred strains.We also noticed that TU_SRA rather than AB_SRA was phylogenetically close to M-AB and *AB, potentially suggesting that AB fish have highly diverged by longterm maintenance in institutions.

Genomic characteristics of the M-AB and IM strains
To address the genomic features of the M-AB and IM strains, we evaluated the genomic homogeneity and heterogeneity of the inbred fish.Homogeneity was defined by the ratio of homozygous nucleotides (number of homozygous nucleotides/number of nucleotides in whole genome of 1.345 GB).The heterogeneity of the M-AB (0.011 ± 0.002%, n = 3) and IM (0.008 ± 0.001%, n = 3) strains was markedly lower than those of the parental *AB (0.086 ± 0.004%, n = 3) and IND (0.197 ± 0.016%, n = 3) strains as well as those of the AB_SRA (0.504%) and TU_SRA (0.196%) strains (Fig. 3A).Since heterogeneities of zebrafish strains were reported based on the total SNP numbers (number of heterozygous nucleotides/number of SNPs found in the analysis) in a previous study 8 , we also calculated the heterogeneity of each strain for the 24,312,666 SNP positions in a total of 14 genomes.The heterogeneity was 0.628 ± 0.129% in the M-AB strain and 0.469 ± 0.042% in the IM strain, which were also markedly lower than those of the *AB (4.770 ± 0.233%), IND (10.937 ± 0.876%), AB_SRA (27.894%) and TU_SRA (10.797%) strains.These data indicate that repetitive sib-pair mating significantly increased the allelic uniformity of the genome in both the M-AB and IM strains.
Although we analyzed the 27th generation of M-AB and the 47th generation of IM, allelically heterozygous nucleotides remained.To explore whether these "heterozygous nucleotides" are distributed evenly throughout the genome or accumulate at specific loci, we counted the number of heterogeneous nucleotides for every 100 kb in each chromosome (Fig. 3B; Table S3).Intriguingly, a prominent hot spot of heterozygous nucleotides was found at the middle of chromosome 3 in two M-AB individuals but not in the other M-AB_2 fish.Since the M-AB_2 fish that carries this region at homozygous survived to adulthood, the heterogeneity in this region would be removed by further sib-pair mating.We also found several genomic loci where heterogeneity was maintained in three M-AB fish (i.e., at the distal end in chromosome 9), three IM fish (i.e., at the distal end in chromosome 23) or all 6 fish (i.e., at the distal end in chromosome 3).These regions might contain closely located harmful genetic variants that are heterozygous for each different chromosomal allele.S3.

SNPs potentially affecting gene products in the M-AB and IM strains
Our WGS of M-AB and IM fish identified a number of SNPs that differed from the zebrafish genome assembly GRCz11.These SNPs were classified into heterozygous and homozygous SNPs in each fish.Generally, repetitive sib-pair mating promotes the accumulation of "homozygous SNPs" that affect the function of the gene products.Such homozygous SNPs represent genomic characteristics of a strain if the SNPs were commonly found in that strain.To evaluate the effects of these homozygous SNPs, we classified the SNPs into 6 groups by location with regard to exon-intron regions and protein coding sequences (Table 2).In both the M-AB and IM strains, 96.9% of homozygous SNPs were found in either intron or intergenic regions, while 1.4% of homozygous SNPs were found in the 5′-UTR or 3′-UTR.The remaining 1.7% were located in the coding region or 2-bp splicing junctions.We then focused on homozygous SNPs that potentially disrupt genes by creating nonsense codons or disrupting initiation codons, termination codons or 2-bp splicing junctions.The overlap of homozygous SNPs that potentially disrupt genes between individuals revealed that most of the gene disruption was commonly found in three fish of each strain (Fig. 4).The list of the disrupted genes provided genetic features of the M-AB and IM strains (Table S4-6).For example, disrupted genes in the IM strain contained several genes relevant for reproduction, such as cep72 16 , cxxc1 17 , cyp1a1 18 , fbxo42 19 , gsto2 20 , hsf5 21 , pex10 22 , psat1 23 , and rad51c 24 .Among them, gene-disrupting SNPs in cyp1a1 (Fig. S5) and psat1 (Fig. S6) generated nonsense codons closer to the C-terminus of the coding sequences and thus might not affect the function of the gene products.However, the other gene-disrupting SNPs potentially affect the function of the gene products (Figs.S7-13).

Morpholino knockdown, CRISPR knockout and Tol2 transgenesis are possible in embryos of the M-AB strain
Since the M-AB strain fish appears to be useful with fertilization and survival rates comparable to those of the parental *AB strain, we finally tested whether M-AB embryos can be genetically manipulated by microinjection  of morpholino oligonucleotides (MOs), CRISPR/Cas9 and foreign genes.MOs against the chordin gene, which is involved in dorsoventral patterning in embryonic development, were injected into 1-cell embryos of M-AB and *AB clutches and subjected to survival and morphological analysis.The survival rate of MO-injected M-AB embryos (87.5 ± 9.6%) at 3 dpf was comparable to that of MO-injected *AB embryos (92.1 ± 6.5%) (Fig. 5A).More than 80% of MO-injected embryos displayed typical chordin-deficient phenotypes, such as expanded blood islands, U-shaped somites and short tail fins with multiple folds 25 in both the M-AB and *AB backgrounds (Fig. 5B).
To attempt genome editing, we next injected CRISPR single-guide RNA (sgRNA) for the sqstm1 gene 26 and Cas9 protein into 1-cell embryos of the M-AB and *AB strains.The survival rates of injected embryos were comparable between the M-AB (74.8 ± 20.6%) and *AB (73.8 ± 26.3%) strains at 2 dpf (Fig. 5C).Our heteroduplex mobility assay suggested that insertions and/or deletions were generated at the target site of the sqstm1 gene in all tested embryos of M-AB and *AB clutches (Fig. S14).
Finally, we injected the ef1α:EGFP plasmid with tol2 mRNA into 1-cell embryos of the M-AB and *AB strains to test transposon-mediated transgenesis 27 .Although the survival rate of injected M-AB embryos (58.2 ± 17.0%) was slightly lower than that of injected *AB embryos (80.6 ± 6.8%) at 12 dpf, the survival rate of M-AB was more than 50% (Fig. 5D).The GFP gene was detected comparably in surviving M-AB embryos (50.9 ± 22.7%) and surviving *AB embryos (37.3 ± 24.9%) by PCR (Fig. 5E).After raising those embryos to the adult stage, we crossed them with wild-type fish and obtained embryos carrying the ef1a:EGFP transgene at 26.0% in the M-AB clutch and at 38.1% in the *AB clutch.Fluorescence detection of GFP was ascertained in 7.1% (9/127) of embryos obtained from the mating of three pairs of injected M-AB fish.These results suggest that the M-AB embryos are applicable to microinjection at the 1-cell stage and thus that M-AB is a useful inbred strain for the analysis of gene function.

Discussion
The present study describes two genetically uniform zebrafish strains generated by consecutive full sib-pair mating for more than 20 generations.Since this criterion is used as a definition of an inbred strain in mice, rats and medakas, we can declare the establishment of IM and M-AB as inbred strains in zebrafish.With subsequent sib-pair mating, inbreeding of the IM strain has recently reached the 50th generation.Although we have provided IM embryos from the 20th-40th generations to the other labs, the IM strain is hardly maintained out of our lab, probably because IM fish lay only a small number of offspring during a short reproductive period.On the other hand, the M-AB strain, which was originated from *AB, did not show such a reproduction insufficiency.We have distributed M-AB embryos at the 26th-28th generations to several labs, and M-AB fish have been successfully maintained over the generation.Since *AB was derived from selected gynogenetic females of the AB strain 12 , genetic variants that adversely affect development and reproduction might have been removed from *AB.In accordance with the elimination of harmful variants, the embryonic lethal TM phenotype, which we often observed in IM offspring and TU-derived sib-pair progeny 11 , was not observed in M-AB embryos.The M-AB strain can be maintained by regular care and thus will be in demand in the zebrafish research community.
Sex is determined by the sex chromosome in mice, rats and medakas, whereas it is determined by both genetic and environmental factors in zebrafish 28 .Despite enormous efforts in sex determination studies, genetic factors that determine sex in zebrafish have not been unveiled.Interestingly, M-AB offspring have a tendency to become male in our breeding.Similar significant male bias has also been observed in the IM strain and in a previous sib-pair mating started from the WIK strain and wild Bangladesh fish 11,29 .Thus, inbred breeding likely promotes male bias in zebrafish.Adjustment of environmental factors enables changes in the sex ratio in zebrafish progeny.In our case, the tendency of male bias in the M-AB offspring was reduced by changing the larval diet from paramecia to rotifers and by decreasing the number of growing larvae and juveniles in a tank.Since the M-AB strain carries possible genetic factors for male bias but produces both females and males, it will also be useful to study genetic and environmental factors in sex determination.
Our previous genomic analysis of the IM strain using 100 SSLP markers suggested that heterogeneity of the IM strain from the 14th generation was 5% 11 .The current study using 24,312,666 SNPs, which we found as a difference from the zebrafish genome assembly, revealed that heterogeneities of the IM strain from the 47th generation and the M-AB strain from the 27th generation were 0.469% and 0.628%, respectively.Apart from these inbred lines, the heterogeneities of AB and TU were 27.894% and 10.797%, respectively.Our study was consistent with a former study using 2120 SNPs, in which the heterogeneities of AB and TU were calculated as 24.8% and 14.6%, respectively 10 .These data lead to the conclusion that the two inbred strains IM and M-AB, which were generated through repetitive sib-pair mating, were genetically highly homozygous.
The zebrafish genome assembly GRCz11 was generated by sequencing the TU genome.Based on this genome assembly, gene information, including transcript isoforms and protein coding sequences, was annotated in the Ensembl database.Our WGS of the IM identified a number of homozygous SNPs, which were shared in three IM fish, showing specific characteristics of the IM strain.The Ensembl suggested that some homozygous SNPs in the IM create nonsense codons or disrupt initiation codons, termination codons or splicing junctions, thereby leading to gene disruption.Such potentially disrupted genes included reproduction-related genes reported in zebrafish or in other vertebrates.For instance, a SNP in the hsf5 gene, which is needed for spermatogenesis in zebrafish 21 , created a nonsense codon at the 357th amino acid residue (Fig. S7).Similarly, SNPs in the cxxc1 and rad51c genes, both of which are essential for meiosis in mice 17,24 , resulted in protein truncation in the middle and first halves of the gene products, respectively (Figs.S8, S9).Other SNPs also suggested the disruption of cep72 16 , fbxo42 19 , gsto2 20 , and pex10 22 , which are relevant for reproduction.Such gene-disrupting SNPs may account for inbreeding depression, such as the low fertility rate of the IM strain.However, we also have to keep in mind that some annotations of genes in Ensembl comprise intron-containing transcripts and thus wrong protein coding sequences due to annotation errors 30 .Some cases of gene disruption retrieved from such homozygous SNPs may be false disruptions.Future studies of loss-of-function of those disrupted genes will elucidate the actual genetic causes of the reproduction-insufficient features of the IM strain.
On the other hand, we could not find homozygous SNPs that affect reproduction-related genes in the M-AB strain.A SNP that disrupts the termination codon of the reproduction-related gene cep128 31 was found in all three M-AB fish and three IM fish (Table S6).However, it is unclear whether this SNP disrupts Cep128 function for the following reasons: the transcript of cep128 (ENSDART00000165489.2) is incomplete, as it lacks a start codon, suggesting the presence of other splicing variants.The M-AB strain is maintained by regular care, and the M-AB embryos are applicable to microinjection for various purposes.Since M-AB is genetically highly uniform, this inbred strain will be useful, especially for the study of behaviors, microbiomes and drug susceptibility, which are often affected by individual genetic differences.Both the IM and M-AB strains are available from the National BioResource Project, which is a zebrafish stock center in Japan (https:// nbrp.jp/ en/ resou rce/ zebra fish-en/).

Animals
The zebrafish *AB strain was kindly provided by Prof. Uwe Strähle in 2010.The IND strain fish was kindly provided by Dr. Yasuyuki Kishimoto.Both the *AB and IND strains were maintained by mass mating of 6 females and 8 males.The heterogeneous fish used in the analyses were bred under the same conditions as the M-AB strain.
Continued inbreeding of the IM strain and establishment of the M-AB strain were performed by basically the same procedure as described previously 11 .To reduce the time for maintenance, some modifications were made after the 28th generation of the IM strain and for the M-AB strain.First, the mating pair that did not produce a large number of offspring was not used for inbreeding.Second, the numbers of eggs and embryos were not exactly counted.Finally, the number of branches for inbreeding was not always five, especially when a sufficient number of offspring was obtained.The best mating pair was determined by subjective impressions of the offspring, including hardness of chorions and morphological normalcy up to 4 dpf.The embryos were raised at less than 50-100 embryos per 150 mm dish with replacement of approximately 1/5 of the system water every day.Feeding rotifers was started at 4 dpf, and brine shrimp were started at ~ 8 dpf, followed by observation of food eating.Five days after starting to eat brine shrimp, they were transferred into the circulation system (~ 20 larvae per 2 L tank).As M.S. has moved to Keio University in 2014, the M-AB strain after the 11th generation has been maintained as a substrain at two places: the National Institute of Genetics and Keio University.
To measure body size, zebrafish were anesthetized in 0.01% ethyl 4-aminobenzoate (w/v), and the body length and standard length were measured using a Vernier caliper.

MO injection
For microinjection of antisense MOs, 2.3 nL of 0.05 mM chordin MOs or standard control MOs (GeneTools, Philomath, OR, USA) were injected into 1-cell stage embryos.Phenotype and survival were evaluated at 2 and 3 dpf, respectively.

Genome editing
CRISPR/Cas9-mediated genome editing was performed based on published protocols 32,33 .Template DNA for sgRNA synthesis was prepared by amplification with a primer specific to sqstm1, a universal reverse primer and T4 DNA polymerase.After purification of the template DNA, sgRNA (taatacgactcactataGGA CGG GTG TGA TGG GCC GGgttttagagctagaa) was transcribed using a MEGAscript T7 Transcription Kit (Ambion, Austin, TX, USA) and purified with a MEGAclear Transcription Clean-up Kit (Ambion).The M-AB and *AB embryos were injected at the 1-cell stage with 2.3 nL of a mixture of 5 pmol/μL Cas9 Nuclease NLS Protein (Applied Biological Materials, Richmond, Canada) and 50 ng/μL sqstm1 sgRNA.We evaluated the mutation at 1 dpf by heteroduplex mobility assay by PCR using the following primers 34 .Forward: 5′-GTC ATA TGG GTC CAT CTC CAAT-3′; Reverse: 5′-AGT GCT ATT CAC CTC AAA CACG-3′.Amplified PCR products were separated by polyacrylamide gel electrophoresis.

Tol2-mediated transgenesis
To generate transgenic lines, microinjection was performed in 1-cell stage embryos with 2.3 nL of injection solution containing 15 ng/μL pTG-egfp 35 and 25 ng/μL tol2 transposase mRNA 27 .Injected embryos were raised to the adult stage and crossed with the wild-type strain to generate F 1 progeny.The expression of EGFP was visually assayed at 7 dpf.Transgenesis was confirmed by PCR using the genome extracted from the caudal fin crip of F0 fish at ~ 30 dpf and the F1 embryo at 1-4 dpf.The following primers were used to amplify the EGFP gene.5′-ACC ACA TGA AGC AGC ACG ACT-3′; 5′-CTT CTG GTT GGG TCT TTG C-3′.

WGS library preparation and sequencing
The genomic DNA was converted to a next-generation sequencing library with 5× WGS Fragmentation mix and Ligase Mix (Enzymatics, Beverly, MA, USA).A total of 1 μL of genomic DNA was mixed with 1 μL of 10× Fragmentation Buffer, 2 μL of WGS Fragmentation Mix and 6 μL of nuclease-free water.The following program was used for fragmentation: 32 °C for 6.5 min followed by 65 °C for 30 min and a hold at 4 °C.Immediately after fragmentation, 4 μL of 5× Ligation Buffer, 2 μL of WGS Ligase and 6 μL of nuclease-free water were added.The following program was used for ligation: 20 °C for 15 min followed by a hold at 4 °C.The libraries were purified twice using 0.8× volume of AMpure XP beads (Beckman Coulter) and eluted in 12

Figure 1 .
Figure 1.Pedigree of the M-AB strain.(A) Pedigree of the M-AB strain from generation 20 to 30.Pairs in italics indicate pairs that are connected to existing pairs.(B) Female and male fish of the M-AB strain.

Figure 2 .
Figure 2. Phylogenetic tree of the M-AB and IM strains.Phylogenetic relationships of the M-AB, *AB, IM and IND strains with open data for AB_SRA and TU_SRA.This tree used 10,000 bootstrap replicates.Coefficients represent bootstrap values for tree nodes.Note that inbred strains have shorter branches than the other strains.The length of the branch indicates the diversity of that strain.

Figure 3 .
Figure 3. Genomic characteristics of the M-AB and IM strains.(A) The percentages of heterozygous nucleotides (number of heterozygous nucleotides/number of nucleotides in whole genome).Error bars represent the mean ± sem.Plots in the bar graph represent biological replicates.(B) The percentage of heterozygous SNPs in every 100 kb region from chromosomes 1 to 25 is indicated for three M-AB and three IM fish.Areas of interest are marked by triangles.Closed triangles: heterozygous region in two M-AB fish but not in the other M-AB individual; open triangles: heterozygous region common to all M-AB and IM fish, gray triangles: heterozygous region common to only M-AB fish or only IM fish.Detailed count data are available in TableS3.

Figure 4 .
Figure 4. Number of genes potentially disrupted in the M-AB and IM strains.The homozygous SNPs that create nonsense codons or disrupt initiation codons, termination codons or 2-bp splicing junctions within protein coding genes were evaluated as gene disrupting SNPs.The upset plots indicate the overlap of disrupted genes in three M-AB and three IM fish.

Figure 5 .
Figure 5. Survival of M-AB embryos was hardly affected by microinjections.(A) Survival of chordin MO-injected M-AB and *AB embryos (n = 5).(B) Approximately 80% of chordin MO-injected embryos showed chordin-deficient phenotypes.The numbers in the color bars indicate the number of samples.(C) Survival of M-AB and *AB embryos injected with sqstm1 sgRNA and/or Cas9 protein (n = 4).(D) Survival of M-AB and *AB embryos/larvae injected with ef1α:EGFP plasmid and tol2 mRNA (n = 9).Plots connected by lines indicate the same experimental group.(E) Percentages of fish (30 dpf) that harbor GFP-transgene detectable by PCR (*AB: n = 3, M-AB: n = 7).Error bars indicate the standard deviation.
μL of 5 mM Tris-HCl (pH 8.0).The adapter for the ligation step was prepared by annealing 100 mM of C*A*C*TCT TTC CCT ACA CGA CGC TCT TCC GA*T*C*T and /5Phos/G*A*T*CGG AAG AGC ACA CGT CTG AAC TCC AGT*C*A*C (*

Table 1 .
Number of SNPs in each fish.

Table 2 .
Number and ratio of homozygous SNPs.